#!/usr/bin/python3
# coding=utf-8

import numpy as np
import os

def logit(x):
    """Logit函数：计算ln(x/(1-x))，带数值安全处理"""
    # 安全裁剪范围 (避免x=0或x=1)
    x_clipped = np.clip(x, 1e-5, 1-1e-5)
    return np.log(x_clipped) - np.log(1 - x_clipped)  # 拆分为两个对数运算提高精度

def gen_golden_data_logit():
    # 配置参数（保持原有结构）
    dtype = np.float16
    input_shape = [8, 2048]
    
    # 生成输入数据（范围设为0.001到0.999，避开边界）
    # 使用beta分布生成更合理的数据分布（α=β=2）
    input_x = np.random.beta(2, 2, size=input_shape).astype(dtype) * 0.998 + 0.001
    
    # 计算Logit函数
    golden = logit(input_x).astype(dtype)
    
    # 创建目录结构
    os.makedirs("./input", exist_ok=True)
    os.makedirs("./output", exist_ok=True)
    
    # 保持tiling参数（兼容性）
    tiling = np.array([16384, 8], dtype=np.uint32)
    tiling.tofile("./input/input_tiling.bin")
    
    # 保存文件
    input_x.tofile("./input/input_x.bin")
    golden.tofile("./output/golden.bin")

if __name__ == "__main__":
    gen_golden_data_logit()